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INTRODUCTION 


Ri LI — meme 


The Double Precision Flostins Point Packese is è convenient 
wy and rowerful tool which adds a whole ranse of new features 
l Lo. the- standard. Arrlesoft. BASIC language., Frosrams wou 
write which reauire more than the normal 9-disit accuracy of 
= Arrlesoft can now use the 2i-disit arithmetic carability of 
| qu DFP reckase. Nothing. is. taken away? all. of 
the standard Arrlesoft features are still available. [FFF 
ot simply augments Arrlesoft in a non-obtrusive manner. 


i Here is a summary of the most excitins features: 


Full 21-Disit Precision in Arprlesoft compatible format 


Sa 


ârithmetic Orerations (t*r-»55/) 


arithmetic Exrressions 


|| 


21-disit Constants end Variables 


INFUT and PRINT statements handle 21-disit values 


A t4 


Automatic loader hides DPFFF below your Aprlesoft rrogram 


* 
X 
x. 
* 
* 
€ Simple Interface with Arplesoft using the "&" Feature 
* 
X Faster efficient Assembly Landuase Codins 

* Compacts only 2048 bytes. 

x 


Works with: Arelesoft ROM card» 
Arprlesoft in RAM» or 
Applesoft in the Lansuase Card 


Works with both Arrle II and Arrle Ii Flus 


»* 


* Frosárams are included which comrute these 
standard mathematical functions in double precision: 
ABS» SAR? LOG: EXP» SIN» COS» ATN 


EER | PERAE 1 ae l 


x Continuins gevelorment and surrort by 
to rrovide more useful features and extensions with 


| future relesses. 


| (Applesoft and Apple II are trademarks of the Apple 
Computer Inc., Cupertino, California.) 
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QUE RE FIWA ACKAGE? 


When Diu You Need More Frecision? 


There ace mans Problems which arise in computing: which 
reauire more than nine sisnificant fisures for accurate 
results. 


For examples Civil Engineers tyrically use numbers to 
express earth coordinates which contain more than nine 


significant fisures. 


We need more than nine didils to express bis money fisures: 
with nine digits we can only express uP to $929999999.99, 


Many sciences use matrix algebra to solve Problems» usins' 
such techniques as matrix inversion and eigenvalues. 
Alihoush data for their analyses mas not reauire more than 
nine digits, the internal manirulations dos otherwise the 
results are erroneous due to the effects of round-off error, 


Problems in statistics: such as correlation or least sauares 
curve fittings may need double precision caPabliitw at 
certain stases in ihe commutation. 


Other Possible Solutions: 


Pe ns ee 


It may be Fossible to Obtain double rrecision when vou need 
it in other wavs. 


One metnod is to use a different computer, Since the Arrle 
II is our favorite computers we do not want to even sussest 
anything along this line. In ans case» the cost of usins a 


different comruter can be overwhelming. 


Another method is to bus a8 srecial hardware device which 
plugs into an Arrle I/0 slot» and dives hardware floating 
roint in both single and double rrecision. Not only is such 
8 board cuite exrensive (if available at all)» but also your 


rrograms could not be written in standard Arrlesoft. 


à third possible solution is to rrosram double rrecision 
calculations usins standard Arelesoft arithmetic orerations. 
This can be done» but it is cuite tedious. It is also vers 
slow in execution time. 


In tne judsment of Guam, the best of all Possible 
solutions is to Provide an assembly lansusse rackase which 


interfaces to standard Arrelesoft through the amrersand 
vector. DFFF is Just such 2 package! 
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HOW TO USE DPFP 


Losdins DFFFP;: . 


lA d 


' The S-C DPFF Fackase is automatically loaded below sour 


Arrelesoft program by the DFFF-LOADER. Simely include the 
followings line in your Frosram as the first line: 
10 FKINT CHRKS$( 4)" BLOAD DFFF-LOADER" $CALL768 


The- DPFF-LOADER loads into memory. from $0300 through $03C7. 
The loader works in several Phases! 


1. Determine whether DFFP Packase already loaded. 
If sos return immediately. 


2. Determine which form of Aer lesoft is being used. 
(ROM or RAM) 


3. Move user's rrosram ur in memory 2048 bytes, 
This makes room for DPFFP. 


4. Adjust all Arrlesoft Pointers to this new Position. 


5. Restrins the Arrlesoft rrosram line address chain 
for the new rosition. 


6. Set ur the "&-vector" at $03F5 to Point to DPFF. 


7. BLOAL the OFFF FackKaser between $0800 and $OFFF 
if ROM or Lansuase System 6rrlesofti» or between 
$3000 and $37FF for RAM árrlesoft. 


8. Perform a CLEAR statements and return to user. 


Once the DFFP Package has been Joadeds it will stay in 
memory until you tyre the FF commands the INT commands or 
re-boot DOS. Of courses: you could destroy it by FONE's or 
hy BLOADins something else over the tor of it. 


You can use. the NEW» LOAL? SAUE commands» and enter or edit 
Programs as sou wish: DFFP will remain intact under your 


Programs and you can call it whenever you wish. 


You will probably want to copy the DPFP-LOADER file and the 
BPFP-1.0-ROM (or -RAàM) files to your own disks. Use the 
followins commands to do that. 


JBLOAD DPFP-LOADER 
JESAVE DFFF-LOADERsA$300rL $C8 


JBLOAD DPFF-1.0-ROM 
JESAVE DPFF-1.0-ROM»6$800;:L$800 
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The DPFF Package is called by the &DF; statement. The 
amrersand is a srecial feature of Arrlesoft BASIC (described 
on rase- 123. of; the àrrlesoft. Reference: Manual ).. This. 
symbol» when executed as an instruction» causes an 
unconditional Jump to location $03F5. The DFFP Loader sets 
up this location to Jumr into the DPFP Fackasee At this 
point the DPFF FackKase bedins to analyze and execute the 
characters which follow the amrersand. 


N 


DPFP will return a syntax error and stor the rProsram if the 
next three characters after the amrersand are not "DF". If 
the characters are correct: DFFP will Proceed to analyze 
succeeding statements serarated by colons on the same line. 


150 &DP:INFUT X(0) 
160 &DP:YCO) = XCOD) X X(OD & 3.14159265358979 323846 


SPRINTY( 0) 


It is not necessary that the "&DFji" be the first statement 
in a line. For examrler the following statement Performs 
the ABS function on a double Precision variable. It uses a 
single Precision comparisons and only the ortional 
complement oreration is in double Precision. 


170 IF X«0) < 0 THEN DF! X(0) = -X(0) 


Two colons in a row with no statement between them will shut 
off DPFF and revert to normal ArrPlesoft processing. 


180 IF X(0) > 1 THEN ADF? X(0) = ACO) X X(O)D 33 GO TO 
100 


You can also ture double Precision statements as direct 
commands in Arrlesoft» once the DFFF Fackase has been 


loaded. 
J& DF? INFUT XCOO0: XCOO02XC 0 KXCO SPRINT X(O) 


(Note that in normal Arrlesoft» the above statement would 
break two rules! multiple statements separated by colons 
are not normally allowed» and the INFUT statement is 
normally not a legal direct statement. Nevertheless» it 
does work in DFFF ausmented Arrlesoft.) 


Study these examrles and the ones in Charters 4 and 5» and 
try your hand at rrosrammins a few yourself. 
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Ail verisbles referenced by DPFF must consist of two 


addecenit arras elements. The array must be a REAL array 


(inal iss the erras must not be INTEGER or STRING>. The 
ArPlesoft Reference Manual describes. arrays. on. Pages. 14». 15» 
aid 98. 


Remember that Arrlesoft acras subscrirts besin with 0 and go 
ur to the limit defined in the DIM statement. An array 
Gimensioned "3r10510" hes three dimensions. The first runs 
from 0 to 35 the second from 0 to 11% and the third from 0 
to li. It could contain 4*12X12-576 real elements» or 
2&l2*122288 double precision elements. 


üárrlesoft arrays sre stored in memory with the leftmost 
eubscriri varsing tne fastest. For examples in the arras 
XYC3»10»10) element XYCO»JsK ) comes immediately before 
element XYCl»JsK 0 Therefore sou mau» in effecti» create an 
array of double precision values by merely Prefixing an 
extra dimension to ihe dimension list. 


If sou wish to set ur serarate variables» you do so by 
dimensioning them to have two real elements each. For 
examr les the statement 


10 DIM ACL )eBC1)xCC1 XC 1) 


will set ur four eerarale variables for use with DPFF. You 
reference the variables within double Precision statements 
with the subscrirt 0. For examrle» 


20 ALFK 0 =AC 0 4 ECOOXCCO) 


If you wish to create an arras of double precision values» 
wou do so bs dimensioning the arras with one extra 
dimension. The extra dimension comes first» and should be 
"i"; this gimension generates two real items» or one double 


Precision item. For exame ler 
10 DIM AC 1si2)rRi1rS26) 


creates bwo arrays which can be used for double Precision 

values. The array A can be thought of as an arras of 13 

double precision values? from AC020) to A Ori2). The array 

B could store 42 double precision values: from RC Os0»s0)2 to 

B.0»5»6)., If you alwass remember to use one exira 

dimension? to rut that extra dimension first? to set that 
"1"5 and to refer to items with the first 


dimension to $ 
subscript = Os then you will succeed in usins DFFF. 
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Constants’ 


Douvle rrecision constants are entered in the same way that 
Single Precision constants are. See pases 4 and 5 of the 
Arrwlesoft Reference Manual for a descrirtion of the legal 


formats. 


The only difference between standard Arrlesoft and the DFFP 
constants is thet DFFF converts and stores 21 sisnificant 
disits» rather than 9. l 


Constant conversion is a time-consuming process: especially 
in double precision. If you are interested in speed» then 
sou should not use constants inside a loor that will be 
executed a lot of times. Instead» store the constant values 
in variables» and refer to the variables inside the loor. 


* 


Conversion between Single and uble Precision: 


It is freauentis desirable to be able to convert a single 
precision value into a double precision one for Purroses of 
computation. This is most easily done bw merely assisnins 
the value to the first cell of a double rrecision variables 
and settings the second cell of the DF variable to zero. 


100 REM CONVERT X TO DOUBLE PRECISION VALUE © 


110 DIM BPC 1) 

120 INFUT X 

130 :$ DFCOOZX? TFC 1 =0 
140 AIF? PRINT DFCO) 

150 GO TO 120 


It is also vers useful to be able to convert from double 
Precision back to single precision. If you only need the 
most sisnificant rart of the double Precision value» then 
there is no conversion to be done? merely reference it and 
you have it. 


However» if sou wish to get at the least sisnificant rFaerti» 
that is also rossible. The followings rrosram shows how; 


100 REM GET SECONDI HALF OF DPFP VALUE 

110 REM AS A SINGLE PRECISION VALUE 

120 TIM X(1)0»T(l1) 

130 PRINT "X = "$i&DPiINFUT X(0) 

140 TCO)0zXCO)0$3 TC 12070 

150 &DFiTCOOZXCOO-TCO) 

iéQ FRINT "LEAST SIGNIFICANT HALF = "sT(CO) 
179 GO TO 130 
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Arithmetic Expressions: 

You can use simple arithmetic exrressions in DPFP. In order 
to Keer the memory reauirements to a minimums some of the 
usual features you exrect in arithmetic exrressions are 
absent. However» enoush rower is present that. you can still 
solve any Problem sou wishe 


Only the four arithmetic operators (+r ~s X» /) can be used. 
Only double Precision variables and constants can he used as 
operands. No function references are permitteds and 
parentheses are not allowed for grouping. Orerations are 
done in a strict left-to-right order. 


The following examples illustrate the carabilities DpPFP 


Provides. | 


10 &DP:XCOOSACODEBCOOXCCO) 

DPFP will first add ACO) and RCO)» and then multiply 
the sum by C(O). In normal Arplesoft exrressionss; the 
precedence rules would evaluate the Product of B(O) and 
C(O) first» and then add ACO) to the result. The two 
orders give entirely different results. 


20 &DP:XCO02-XCOD 
Unars minus is allowed. 


30 GDP STC O =XC 0412 YC O)7XC O)-1/TC O) 

Computes y = (x-1) / (x41). Since DPFF does not 
support parentheses in expressions, you sometimes have 
to resort to using Lemrorars variables such as T(O) and 
compute the expressions in ranis. 


40 &DP:PCOD) = KCOSS) X XCOD + KC0s4) X XCO) + KCO»3) X 
XO) + KC Or2) K XCO) + KCOs1) X XCO) + KC O20) 

Computes the fifth-desree Polynomial rPzkKSXxf9 t K4Xxf14 
+ k3%xł3 + K2Xxf2 t KlXx + KO 


You must remember that DPFP does not have any effect on the 
interrretation of ordinary sinsle Precision expressions. It 
only affects expressions after, and on the same line with, the 


"&DP:" which calls DPFP. 
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INPUT Statement: 


The DPFP Package surrorts a limited fors of the INPUT 
statement, The limitations are these? only double 
Precision variables can be used in the input list» and no 
prompting string im euoles can: be used. 


One double Precision variable» or a list of double precision 
variables sePareted by commas can follow the word INPUT. 


10 &DPS INPUT ACO) 

20 &DP:INPUT ACO)» BCOsN) 
30 FOR I=0 TO 4 

40 &DP:INPUT BC OT) 

50 NEXT I 


If you wish to prompt the user» wou can do so with a PRINT 
Statement, 


10 PRINT "ENTER R("$;IS")"$ 
20 &DP:INPUT X(O»I) 


When the INPUT statement is execuled: a Question mark will 
be displayed and the Aprle will wait for sou to enter a 
numeric value. The format for input numbers is the same as 
for standard Applesoft. Of courses you sav ture ur to 38 
digits» Just as in standard Arrlesoft. The difference is 
that standard Arplesoft will only store 9 significant 
disits» whereas DPFP will store 21 sisnificant digits. 
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PRINT Statement: 


A d 


The DPFP Package suprorts a limited form of the PRINT 
statement. The limitations sare these: Only double 
precision variables or constants mau be in the list of items 
to be Printed? only commas meu be used to serarste the items 
im the print. liste. i 


Since each double precision value may take uP to 27 
characters to Prini» onis one value will be Printed Per 


line. 


All values. are Printed in the general format: 
sd. dddddddddddddddddddGEs dd 


The first "s" is nothing for Positive numhers» or "=" for 
negative numbers. The first "d" js the first significant 
didit (never 0 unless the whole number is zero). Ur to 20 
more digits will be rrinted. Trailing zeroes will not be 
printed. If the value has a decimal exponent other than 
zeros the letter "E"» a + or - signe and s two digit 
exponent will be printed. 


0 


1.5 
6«66666666666666666666E-01 
3.14159265358979 323846 


The format is almost identical to that used in standard 
Applesoft. The difference is that standard Arrlesoft will 
Print values which contain the decimal roint answhere within 
the nine significant fisures without a decimal exronent. 
DPFP will always use the scientific notation form. 
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PEEKing and PONEins in Variables: 


A special feature of DPFP is the ability to obtain the 
memory address of the first bste of a variable. This 
ability allows you to Perform special functions not Possible 
with standard mathematical orerations. 


The feature is invoked by using the "5" symbol as an 
orerator to store the address of a double precision variable 


into an inteser variable: 
10 AZ=O3&DP 3 XC 0 -AZ 


will place the address of the first byte of the value of 
X(0) into the the integer variable AZ. Why was the variable 
AZ first set to 0? Because if the variable AZ had not wet 
been used» Arrlesoft would have created it after fisurins 
the address of X(0). This Process of creating a new inteser 
variable changes the address of the double precision 


variable! 


The address stored in the inteser variable mas be used in 
PEEK and POKE statements. However» it is imrortant to 
remember that the address is only valid until you create 
another non-array variable. 


This address location feature is used in the lodarithm 
program given in charter 4. Since the losarithm (hase 2) of 
a number is ecual to the losarithm of the fractional sent + 
the Power of two exrressed by the exronent party the 
logarithm subroutine extracts the exronent part and then 
sets the exronent to 210. 


Another use of the address location feature is shown in the 
following examples which will display the hexadecimal value 
of any double precision number which is entered. 


O PRINT CHR$(4)"BLOAD IPFP-LOADER"SCALL 768 
5 PRINT "DISPLAY DPFFP VALUES IN HEXADECIMAL” 
10 DIM X(1)0:D1-20:D2-20:BYTE-O0:I-O 

20 AXsOI&IP:XCO) > AZ 

30 PRINT "X = “si&DF$ INPUT X(O) 

40 FOR I = 0 T0 9 

50 BYTE = FEENCAXTIO:GOSUB 500 

60 NEXT I 

70 PRINT 

80 GO TO 20 


500 Di = INTC BYTE/16)3 D2 = BYTE - D1X*16 
510 IF D1»9 THEN D1=fi+7 

520 IF f2>9 THEN D2-D247 

530 PRINT CHR$CD1f4480$ CHR$CD2t480" "5 
$40 RETURN 
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STANDARD MATH FUNCTIONS © 


a ERREUR 


The standard mathematical functions ARS: SRR» LOG» EXP» SIN» 
COS» and ATN are not included in: the: DPFP PacKase. However; 
they can be easily programmed in ArPlesofi: using the DPFF 
carabilities. 


Programs for each of these functions are included on the 
DPFF Diskette. Thes are also listed in this charter» for 
«our study as examrles. You may easily include them as 
subroutines in your own larger Applesoft prosram. 


Alsorithms for the LOG» EXP» SIN» COS» and ATN functions 
were taken from the book COMPUTER APPROXIMATIONS: by J. F. 
Hart et al» Published by John Wiley & Sons» New York» 1968. 
LOG» SIN» COS» and ATN are arrroximated with the same 
general techniaue: 

1) check for valid argument? 

2) reduce the range of the arsument; 

3) using the reduced argument z» arr-roximate F(z) = 2 X 

P / Qs where P and Q are rolynomials in z125; 

4) manipulate F(z) appropriately to expand the range 

back to F(X). 
EXP is approximated in similar stersr excerpt in ster 3; 
F(z) = (zkF*Q0/CzXP-0). 


ARS -- Absolute Value’ 


The absolute value function is simple enough that you can 
code it in a single line: 


100 IF X(0) < O THEN &DP2X( 0 )=-x(0) 


SQR -- Sauare Root: 


i — PS ———— cmmm—À—— 


The techniaue used for the SRR function:,is to first use 
standard Arrlesoft SQR for a sinsle Precision arrroximatian 
to the double Precision seuare root. Then use the 
Neuton-Rarhson iteration twice to improve the Precision to a 


full 72 bits. 


O PRINT CHR$(4)"BLOAD DFFP-LOAIIER" 2CALL768 
5 REM SQUARE ROOT 

10 DIM X(125»8C1) 

20 PRINT "X = "$3&DP:INPUT X(O) 

30 GOSUB 100 

40 PRINT "SQRT(XOzZ"$i&IP;PRINT SCO) 
50 GO TO 20 

100 SCO)OzSQRCXCO 2218612070 

110 &DP:8CO)02XCO/8C 0 481 0 )/2 

120 &DF:8S(0)02XCO02/8CO 4800/2 

130 RETURN 
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LOG ~- Logarithm: 
The logarithm function is computed by the Prosram belows 
using alsorithm #2706 in HART. 


To include the losarithm as a subroutine in your Arr lesoft 
programs you need lines 10-80 to set ur the constants for 
the algorithm, Lines 500-590 are the subroutine. You will 
probably want to chanse the error return in line 500. 


Q PRINT CHRS (4)"BLOAD DPFF-LOADER":? CALL 768 


9 REF 
DIM X(1)0»8C10» TCL) KFZ )» KGC 9 )s PC 22s QC 1 2» ZC 1 DV EC 3)»K2C 3) 


COMPUTE LOG(XO» #2706 IN HART 
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10 

20 à DP:kKFCO) = 5.13904588649239920696E + 02 
25 & DPikF(2) = - 7.92102505773443199064E + 02 
30 å DPSKPC4) = 2.40707633649031186636E + O2 
35 & DP:KF(6)- - 3.5419160305337449948E + O1 
40 å DPSKQCO) = 178.1057583495195620388929 

á5 & DPIKQ(2) = - 333.8903954121714992817202 
50 & DP:KQC4) = 193.7559146303587951736132 

55 & DP:NQ(6) = - 325.5262511104002387 351584 
60 & DPSKACB) = 1 

70 & DP:K2CO) = 1.4142135623730950488 

80 å DP:K2C(2) = .6931471805599453094172321 


110 PRINT "X = "$: & DP? INPUT X(0) 

120 GOSUB 500 

130 PRINT "LOG(X) = "$i & DPS PRINT ECO) 
140 GOTO 110 


S00 IF X00 < = 0 THEN PRINT "RANGE ERROR IN LOG" 
FOF i GOTO 110 
$10 AX = Of & DPi2(00 = XCODiZCO) > AZii REM GET ADDR 


ESS OF Z(0) 

520 E(0) = FEEK (AX) - 128:EC1) = Oi REM GET POWER 
OF TWO 

530 POKE AX»128: REM MAKE Z(0) A FRACTION 

540 & DFiZ(OOD = ZCO) * K200) 

550 & DP:T(0) = 260) t 132(00 = 2000 - 1 7 TC05:8C0) 
= Z(0) X Z(0) 

260 & DPiP(O) = KPC6) X SCO) + KPC4) X SCO) + NFC2) X 
S(O) + KFCOD X Z(O) 

570 & DP:G0(0) = KACB) X SCO) + KRACK) X SCO) + KQC4) x 
S(O) + K2) * S(O) + KACO) 

580 .& DP:ECO) = FCO) / QCOD + ECO) - S X K2¢ 2) 


990 RETURN 


ap = 


r 
Festum Mp cale etae V =e a maa md P 


MU ER OET 


a CNW Ims 


Lol 


EPUM PATI eer rrt T ese, Sce n me n CREARI te 


! 





EXP -- Exronential: 


The exponential function (etx) is computed by the rrodram 
belowr using alsorithm #1069 in HART. 


To include the exronential as a subroutine in your Arrlesoft 
prosremr you need dines i0-95 Lo: set ur- the» constants for 
the algorithm. Lines 1000-1110 are the subroutine. You 
will probably. want to chanse the error return in line 1010. 


Q FRINT CHR$ (4)"BLOAD IPFF-LOADER": CALL 768 
S PRINT "COMPUTE EXFOXO» USING #1069 IN HART” 


20 
30 
40 
90 
60 
79 
80 
90 
95 
100 
110 
120 
130 


1000 
1010 


1020 
1030 
1040 
1050 
1060 
1070 


1080 
1090 


1100 
1110 


pr pr pr or p» pe p p» Or 


IFIKFĖ( 0) = 
TFS KF< 2) 

DFiRFCÁ) 

DPF:RQCO) 

PIKA?) 

DFS KAC 4) 

DPS KRG ó) 

TFS L200) 

TP SL 2¢ 2) 

FRINT "X = 
GOSUB 1000 
PRINT "EXF( 
GOTO 100 


58 " » u u HW p n 


& CFPO) 
IF Z0) > 
FOF’ 
TO) = INT 
T2) = 2 ft 
& DP:2(0) 
& DP:2CO) 
& DP:8CO) 
& DPiF(O) 
x Z(0) 
& DP:0(O) 
* S(O) t 
& DP:T(O) 
& DP:ECO) 
RETURN 


10 DIM X(10»8CA de TC S Ds RECS Os RRC)» PCi )*Q0C 10:ZC 1)» L2C 3) 


2080384.346694663001443B4 
20286.97169744036299076049 
60.61485330061080841615585 
6002720.360238832528230908 
327725 .1518082914423057964 
1749.287689093076403844945 


1 
.6931471805599453094172321 
1. 41421356237 30950488 
"33 & DP; INPUT X(0) 
"$i PRINT ECO) 


X) = & DFS 


= X(0) / L2C0)33 REM CHANGE TO 21Z 
128 THEN  FRINT "EXFONENT TOO LARGE": 


t GOTO 100 - 


(ZC009:T(12 = O 
T(OXTC32 = 0 
= Z(0) - T(0):} REM FRACTIONAL PART 
= Z(0) - .5:1 REN NORMALIZE 
= Z0) & Z(0) 
= KP(4) X SCO) + KFC2) X SCO) + KPCOD 


K0C6) X SCO) + KQC4) X SCO) + KQC2) 
KAC O) 

QCO) - PCO) 

KO) + POO) / TCO) k L2C2) X TC2) 


m 
= 
-m 
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SIN COS -- Sine and Cosine: 


The SIN and COS functions are computed by the program he lows 
using algorithm #3372 in HART. 


To include them in sour Applesoft Prosram» you need lines 
10-80. to set ur- the constants. for the alsorithm. Lines 
500-600 are the subroutine. Notice that the subroutine 
really is for the SIN function only. The techniaue for 
generating a COS is to add Pi/2 to the argument before 
calling the subroutine. (See line 130.) 


0 PRINT CHRS (4)"BLOAD DPFP-LOADER": CALL 768 
5 FRINT “COMPUTE SIN (XO» #3372 IN HART" 
i 


O DIM X(12»8610» TC)» KPC 11 2s KQC9 2sPC 2250€ 12» 2C 1 29 P2 
(3) 

20 DF SKPCO) = 50353932.76969750652546526 

2 DPiKFC2) = - 18889439.03047503439828955 

30 DFiNPC4) = 1838302.585787929312472698 

3S DPiKPC6) = = 70370.60357992078102737647 

40 DPiNPC8) = 1181.360759554700055721965 


DP:NFC100 = - 7.546627236269 5250200665401 


^ 
et 
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50 DP:NQCOD = 22056309.21765732134983309 

55 DP:KQC2) = 1197238, 397665401356441 6306 

60 DP:KQC4) = 19853,25236390720767652221J 

65 DEF:NQCA) = 197.5688875569672025108402 

70 DP;NGCB) = 1 

80 DPSP200) = 6.283185307179586476925:F2( 2) = F2CO 
)/ 4 

100 PRINT "X = "s? & DP: INPUT X(0) 

110 & DPF:Z(0) = X(0)0:3 GOSUR 300 

120 PRINT "SINCX) = "$i & DP; PRINT T(O) 

130 & DP:2(0)0 = X(0) + P2¢2)33 GOSUB 300 

140 FRINT "COS(X) = "$i & DPS PRINT TCO): 

150 GOTO 100 

500 SF = 03 IF Z(0) < O THEN SF = 12 & DP3Z(0) = -Z 


C0) 
510 & DF:2Z(0)0 = Z(0) / Facog) 
520 TCO) =. INT (ZC0091T(10 = OF å DP:2(00 = 760) - T 


(00 X 4 

530 IF 20) > = 2 THEN SF = NOT SF: & DP?Z(O0) = 2 
0)-2 

540 IF Z(02)? = 1 THEN & DPSZ(G) = 2 - Z(CO0) 


$550 å DF:S(0O) = Z(OD kK Z(0) 

560 A DP:8Q(0) = KQC8) X SCO) t NQCAD) & SCO) + NQC AD X 
S(O) + KQC2) * 8COD + KACO) 

570 & DF:iFCO) = KPC10) X SCO) + NPCB) X SCO) + NPCÓD 
X SCO) + NKFC4) & SCO) + NFC2) X SCO) + NPCODO X 


Z(0) ; 
$580 & DFIT(0) = FCO) / QO) 
590 IF SF THEN & DFIT(0) = - T(O) 
600 RETURN 
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ATN -- Arctansent: 


The ATN function is computed by the program belows using 
algorithm #5102 in HART. 


To include the ATN function in your Arelesoft programs vou 
need lines. 10-190: to. set ur- the- constants for the. alsorithm. 
tines 500-590 are the subroutinee The subroutine returns 
values for the arctangent between -Fi/2 and Fi/2. 


O PRINT CHR$ (4)"BLOAD DPFF-LOADER": CALL 768 
© FRINT “COMPUTE ATN (X)» #5102 IN HART" 


10 DIM X( 10» TC10» KPC 17 De KQC 17 De PCL QC 22»ZC 1) y P2C 1)» 
SC1) 

30 & DPSKPCO) = 21514.0596260244193319 3254 

40 & DPSKPC 2) = 73597,43380288444240814981 

50 & DP:KFC4) = 100272.5618306302784970512 

60 & DPiNPCA) = 69439,29720032252337059766 

70 & DF:iNP(8) = 25858.09739719099025716568 

80 & DPiNF(100 s 5038.6391855012666579 37791 

85 & DP:NFCI20 = 460.1588804635351471161727 

90 & DPSKPC14) = 15.08767735870030987717456 

95 & DP&KP( 16) = .07523052818757628444510730 

100 & DP:KQ0(0) = 21514.059626024419331935298 

110 & DP&KQ(2) = 80768.78701155924885176713 

120 & DF:KQ(4) = 122892.6789092784776298743 

130 & DP:KQ(6) = 97323.2034905 3555680260424 

140 & DP:NO(8) = 42868.5765204640809 3184007 

150 å IP:kNGQC10)0 = 10401.13491566890057005104 

160 & DP:N0(12) = 1289.750569116110971411460 

170 & DF:KQC14) = 68.51937831018968013114024 

180 & DFiKQ(160 = 1 

190 & DP:P2(0) = 1.57079632679489661 92313217 

200 PRINT "X = "$i & DF: INPUT X(0) 

210 GOSURB 500 

220 PRINT "ATN(OX) = "53 & DP: PRINT T(0) 

230 GOTO 200 

500 & DP2Z(0) = X(O) 

$10 SF = 0: IF X0) < 0 THEN SF = 1: & DP3Z(0) = -Z 
(0) 

520 FF = 0: IF 260) > = 1 THEN FF = 1% & DP?2(0) = 1 
/ Z(0) 

530 & DP:S(0) = ZCO) X* ZCO) 

540 å DP:QC0) = KQC16) X S(O) + KQC14) X SCO) + KACI 
2) X SCO) + KQC 10) X SCO) + KACB) * SCO) + NRCÓ 
) X SCO) + KQC4) X SCO) + KOC 2) X SCO) + KACO) 

550 & DP:P(0) = KPC16) X SCO) t KFC14) X SCO) + KP(1 
2) k SCO) + KPC10) * SCO) + KP(B) X SCO) + KPC6 
) X SCO) + RPC4) X SCO) + KPC2) X SCO) + NPCOD 

560 å DP3T(O) = ZCO) X PCO) / ACO) 

570 IF FF THEN & DPiT(0) = P2(0) - TCOD 

580 IF SF THEN & DPiT(0) = -~ T(OD 

990 RETURN 
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FURTHER EXANFLES 





Comrulgticu ut FI: 


There ore mots Gifferenl ways LO compute en erFrOXimsie 
value of PI. In fects the value is already documented to 
mang wore: dizils tiau uou. or I will ever need. However» it 
mgkes of interesting exercise to write a rrosram which 
generates a soou value. 


It may be interesting to Know that Edward Shanks» in 1873» 
calculated FI to 707 disits by hend! (Unfortunately he 
made an error in the 528th rlace.) In 1955» FI was Known to 
3089 decimal places} in 1961» a computer was used to stretch 
List record aul Lo 1007060 digits. I would not be surrrised 
to learn thet someone has seneraled several million digits 


of PI bs now! 


One of wans rossible infinite series whose sum is PI can be 
derived from the arctangent series: 

arctan(xJ#x- xt3/3 + xt9/5 - xt7/7 + s.: 
The series in this form converses very slowly, but by 
combining it with the formula for the tangent of the sum of 
two angles» we can get an efficient series. 


The following formula was develored by Machin» in 1706. It 
is the formula used bs Shanks in 1873. 


PI/4^ = 4XalnC1/9) - stnC 1/239), 


4XC1/5 - C1/3X A/S IS + C1/5X 17/3015 ~ D 
-( (1/239) - C1/3X 172391 35. ł ...) 


The following rrozram implements Machin’s method. It 
converses so rapidly that in only 15 iterations we senerate 
a full 21 decimal places! 
d 

O FRINT CHR$ (4)0"BLOaD DPFF-LOADER": CALL 768 

S FRINT “CALCULATE PI" 

10 DIM XC 30»sYCSOsFICIOsMCIOSIC1) 

20 á DPiF4(002 0X10) = 1 7 51Y(02 = 1 / 239 

30 & BFiX€O2) = XCO) X XCO01YC22 = YCOO X Y(O) 

40 SIGN = 1:!D(0)0 = iKi) = 0 

5S0 à DFiM(0024/ IO) 

60 å DPiFICO) = XCO) & 4 - YCO) $ hCOD + FICO) 

70 á DF: PRINT FICO) 


80 SIGN = ~ SIGN 
90 IKO) = € ABS (DCO)) + 2) X SIGN 


100 & DFEXCO) = X(0) X X(2) 
110 à DFiYCO) = YCO) X Y(2) 
120 GOTO 50 
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The coefficient of correlation of a set of data havins X and 
Y values is given by the formula: 


n2xu - 2x24 
p 3 e—————— ÁO 


ie. as FR Ee a UC ORE P 
V(nZxf2 — (Zx2012) (nZut2 - (302) 


Altnoush for most data single Precision calculations ere 
sufficient for fisuring the correlation coefficient and 
other statistical rrorerties of datar when the number of 
data items is very larse significant errors can arise. For 
the PurPoses of illustration the following Program is sivens 
it will calculate the correlation coefficient for X and Y 


vectors. 


O PRINT CHR$ (4)"BLOAD DPFP-LOADER": CALL 768 
5 PRINT "CORRELATION COEFFICIENT" 

10 DIM X(1»200»YC1»s20 )5SC 1» 50» TCS) 

20 REM GENERATE THE DATA 

30 FOR I = 1 TO 20 

40 X(Os»I) = Iix(isI) = 0 

£0 Y(OrI) = I + RNE (1) X 21YCA1»ID) = 0 

60 NEXT I 

100 REM CALCULATE THE SUMS 

110 FOR I = i TO 5: & DF:8COSID = Oi: NEXT I 
120 FOR I = 1 TO 20 

130 å DFíSCOs1) = XC OrI) + S(O»1)33 REM SUMCX) 
140 & IP:8€0»2) = YCOrT) + S(0r2)?: REM SUM Y) 
150 & LDPiSCO3) = XCOsIO K X(OsID + S8(053)0:: REM SUM 


(X12) 

160 & DP:S(0s4) = YCOSIO X Y(O»IO t S(0s4):: REM SUM 
CY12) 

170 & DP:S(0»5) = X(O»IO X YCOsIO 4 SCO, 52:1 REM SUM 
( XXY) 


180 NEXT I 

200 REM CALCULATE THE CORRELATION 

210 & DP:T(0) = S(Os1) & SCOv2) 

220 & DP:T(0) 20 X SCO»S) - TCO) 

230 & DP:T(2) $602) * SC Or2) 

240 &.LDP:T(2) 20 X S(Or4) - T(2) 

250 & DP$T(4) $0510 X* 9SCOs1) 

260 & DPiT(2) 20 X S(Os3) - TC 4) X TCS) 
270 T(4) = SAR (T(2)) 

280 & DPiT(4) T2) / 14) + T(42 / 2 
290 & DP$T(4) T2) / 114) + T(42 / 2 
300 & DP$T(O) TCO} / TC 4) 

310 & DP? PRINT T(0) 
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he trix tniultirlication: 
Since many Froblems arise in matrix alsebra reauiring hisher 
precision? the followings Prosram for matrix multirlication 
is presented as an example. You can study it to learn how 
to use UFFP in other Situations requiring. matrix 


computations. 


O PRINT  CHR$ (&)"BLOAGD DFFF-LOADER"? CALL 768 

& PRINT “MATRIX HULTIFLY" 

10 INFUT IsJsh 

20 DIM &C1sI»sJOsBCAs Je KOsCC As I RO 

30 FOR II = 0 TO I! FOR JJ = 0 TO J 

40 PRINT II",;"JJ"; "5i & DF? INFUT ACO:IIsJJ) 

50 NEXT JJ:II 

&0 FOR JJ = O TO Ji FOR KK = O TO K 

70 PRINT JJ"»"NK"j "si & DP: INPUT BCOsJJs NR) 

80 NEXT KK»JJ 

100 FOR II = O TO I 

110 FOR KK = O TO K 

120 & DP:CCOsIISsNK) 

130 FOR JJ = 0 TO J 

140 & DPCCO»IIsNRO = ACOsII JJ) X BCOsJJs KK + C(O» 
II»KRK) 

150 NEXT JJ 

160. PRINT II"s"RK"j "$1: & DP? PRINT C( Os II» RAO 

170 NEXT KK 

180 NEXT II 
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TIMING 





It is difficult to determine exactly how long a particular 
operations such as addition or subtraction» takes, The 
execution time derends on the vaiues of the orerands» soa 
even counting of machine cycles would not sive an accurate 
answer, Furthermore» an isolated timing is difficult to 
relate to timing in actual use embedded within a KASIC 
Prosram. 


The technioue I used to arrive at some fisures for the 
execution time was io construct a loor and execute it 10000 
times. I first timed it with an emrtw loors and then added 
the feature I wanted to time. The difference in the two 
times should indicate the time for the added oreratian. 


Here is the empty loor? 
10 ACO) = 2.718281828459045235356 
20 AC2) = 3$.14159265358979 323846 
30 FOR I = 1 TO 10000 
50 NEXT I 


The empty loor rrosram took 13 seconds to execute. Nividins 
bu 10000» we discover that the averase time for s sinsle 
ster through the loor was 0013 secondsy or 1.3 


milliseconds. 


Next I added this statement; 
40 REM 


The new time wes 16 seconds» for an averade loorp-ster of 1.6 
milliseconds. Do you realize it takes 300 microseconds to 


execute a REM statement? 


Next I rerlaced the REM statement with a normal single 


precision addition: 
40 AC4) = ACO) + C2) 
The execution time was 101.6 seconds. Subtracting the empty . 


loop time leaves 88.6 seconds» or 8.86 milliseconds rer 
execution. 


Next I replaced line 40 with the following DPFP statement: 
40 &DPtac4) = ACO) EAC?) 

The new time was 109 seconds» only 7.4 seconds lonser than 

the single precision version. We det 9.6 milliseconds for 


the DPFF statement. 


Performind similar tests for the multirls and divide 
operations» I comruted the following times: 


AC4) = ACO) X AC?) 10.82 milliseconds 
ale SAC 4 =AC 0 KA 2) 270.08 milliseconds 
ac4) = ACO) / AC2) 11.08 milliseconds 
&DP2406O 4 =AC0)/A0 2) 22,74 milliseconds 
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If we make the assumption that these are typical times for 
any random pair of arguments, we see that wmultielys and 
divide take roughly twice as long in DPFF as in normal 
Applesoft. Addition and subtraction are very nearly the 
same in both precisions. 


The following chart grarhicaliy illustrates the timing 
information. discussed. above., 
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Internal Format: 


Values are stored in ten consecutive bytes in memorwv. They 
occupy two. adjacent cells of an Arrlesoft array. The format 
of the first five bytes is identical to the standard sinsle 
precision floating point format in Arrlesoft. The remaining 
five bytes are merely added precision for the mantissa, 


Floating Point values consists of two Parts? the exponents 
and the mantissa. The exronent in DPFF values is stored in 
the first byte. It is an excess-128 value. This means that 
an exponent of O will be stored as $80 (04128)s 1 is $81» 
and -1 is $7F. The exponent indicates what rower of tuo to 
multiply the mantissa by to get the number's true velue. 


The mantissa» also called the "fraction", is 3 9-hute 
fractional value. That is» the radix Point is assumed to he 
to the left of the first bit of the mantissa. The mantissa 
is kert in normalized position. This means that it is 
shifted to the left until the most significant bit is a 1. 
After normalizins» the leftmost bit is no lonser needed 
(because it is always =1). Therefore we can use its 
Position for the sign bit of the value. 


The first byte is the biased exronent; the first hit of the 


next bute is the sign? the remaining 71 hits» Plus the 
"hidden-1" bit in the sisn bit Position sre the mantissa or 


fraction. 


Mere are some examples: 


80 00 00 00 00 00 00 00 00 00 = 0.5 

80 80 00 00 00 00 00 00 00 00 = -0.35 
81 40 00 00 00 00 00 00 00 00 = 1.5 

00 00 00 00 00 00 00 00 00 00 = 0 

81 OC CC CC CC CC CC CC CC CI! = 1.1 


During comrutation inside ItPFFs an extra bete is carried 
along, The extra byte sijves Protection against truncation 
error. After each orerationr the extra bute is used to 
determine whether the value in the floating roint 
accumulator should be rounded or not. 


293 EA 
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Precision? 


Precision is a function of the internal formal: and derends 
most heavily on the number of bits in the mantissa. You can 
comrute the number of decimal places which a given number of 
bits can rerresent by multiplying the number of bits Dv ihe 
los base 10 of 2. in other words» with 72 bits in the 
mantissa we can represent 72X.30103 = 21.67416 decimal 


disi ts. 
What on earth is ".67416" of a decimal digit? The number 
21.67416 means that we have a full 21 disits of Precision: 


plus a few extra guard bits. The larsest inteser which can 
be exactly rerresented in 72 bits is 2172 - 1» or 


4972223669482 9869 564552133695 


If the number above were stored in DPFF» it would look like 
this in hexadecimal: 


C8 7F FF FF FF FF FF FF FF FF 


Rande: 


Rande is determined by the number of bits allotted to the 
exponent field. Since DPFF allows eight bits for the 

exponents it can indicate any value from 21-127 to 211277. 
This translates roughly to a ranse from 101-38 to 10f438. 


(127%.30103 = 38.23081) 


The smallest positive number which can be represented in 
DPFP is? 


01000000000000000000 = .Sk2t-127 
= 2.93873587705571876992E-39 


Any smaller value which results dur ins commutation will be 
stored as zero. uà 


The largest rositive number which can be rerresented in DPFF 
is: 


,9999,,, 921127 


FFZFFFFFFFFFFFFFFFFF = 
= 1.70141183460469231731E £28 
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v REGISTRATION FORN! DOUBLE PRECISION FLOATING POINT PACKAGE 
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\ 
Please fill out and mail this registration form to S-C 


SOFTWARE at the address above. You will be informed of 
improvements and new versions of the DPFP Fackase as they 


become available. 
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